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Electrochemical etching of silicon in hydrofluoride containing electrolytes leads to pore formation 
for low and to electropolishing for high applied current. The transition between pore formation and 
polishing is accompanied by a change of the valence of the electrochemical dissolution reaction. 
The local etching rate at the interface between the semiconductor and the electrolyte is determined 
by the local current density. We model the transport of reactants and reaction products and thus 
the current density in both, the semiconductor and the electrolyte. Basic features of the chemical 
reaction at the interface are summarized in law of mass action type boundary conditions for the 
transport equations at the interface. We investigate the linear stability of a planar and flat interface. 
Upon increasing the current density the stability flips either through a change of the valence of the 
dissolution reaction or by a nonlinear boundary conditions at the interface. 
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I. INTRODUCTION 

Porous silicon was discovered in the fifties trying to 
electropolish silicon in hydrofiuoric acid For low 

current densities, respectively high electrolyte concentra- 
tions, silicon is not electropolished but pores are formed. 
Increasing the current density over a threshold value, 



which decreases with the electrolyte concentration, re- 
sults in electropolishing. In the beginning of the nineties 
visible luminescence at room temperature was discovered 
The possibihty to produce optoelectronic devices 
out of porous silicon started enormous research activity. 
Meanwhile many applications for porous silicon are in de- 
velopment. Most of these applications are based on the 
morphology of porous silicon, for a review see ||]. 

Porous silicon is formed by anodic dissolution of sili- 
con in hydrofiuoric acid. The silicon surface is in contact 
with the electrolyte, usually in a Teflon cell. Through 
the external potential an electric current is maintained 
across the cell and flows from the semiconductor to the 
acid. Defect electrons (i.e. holes) from the semiconduc- 
tor and HF (or F~ ions) from the electrolyte combine 
at the fluid-semiconductor interface and dissolve silicon 
through an electrochemical reaction. The morphology of 
the unsolved silicon depends on the current. In the elec- 
tropolishing phase the silicon surface is etched layer by 
layer and remains essentially flat, whereas in the porous 
silicon phase many holes are formed of a size ranging 
from a few nanometers to microns. Porosities of over 
95% relative to crystalline silicon can be reached. 

Despite its importance, there is little theoretical under- 
standing of how porous silicon is formed, for a review see 
1^. Even the basic issue why there is a transition from 
porous silicon formation to electropolishing is unresolved. 
The reasons are rather obvious, when one compares with 
other pattern formation processes like dendrites, viscous 
flngering, and colloidal aggregation. As in these systems, 
we have to understand the dynamics of a moving inter- 
face, here between a semiconductor and an acid. At the 



*email: rauscher@ccmr. Cornell. edu 



1 



II MODELING ELECTROCHEMICAL ETCHING 



2 



interface silicon is dissolved through an electrochemical 
reaction. Thus in contrast to the better understood sys- 
tems mentioned above, we have to consider the trans- 
port of several species in the presence of an electric field. 
Moreover, since the species feed back into the electrical 
field and react with each other, the transport equations 
become nonlinear. 

One approach to model the dissolution process is by 
a stochastic growth model, i.e. growth of the fluid into 
the semiconductor, see for example 0,^. These models 
are inspired by the diffusion limited aggregation (DLA) 
model. While, in principle, such a model is on the atomic 
scale, because of numerical limitations in practice larger 
spatial units are used. By suitable adjustments of model 
parameters structures qualitatively similar to porous sil- 
icon can be produced, but the quantitative connection of 
model parameters to physical parameters is lost and in 
some cases the coarsening introduces a new length scale 
into the system which obscures the physics. 

A second approach, starting at much larger length 
scales, is to use a continuum description for both, the 
motion of the interface and the ionic and electronic trans- 
port, see e.g., |9|-|l^]. These models include a depletion 
or passivation layer at the interface in a phenomenolog- 
ical way and assume an ad hoc surface tension to stabi- 
lize against small perturbations and to provide a length 
scale which can be compared with the pore formation. 
However, surface tension can only affect length scales of 
the order of the micropore diameter, i.e. nanometers. 
The above mentioned models are not valid on these small 
scales since the mean free path, quantum effects, and the 
electrical double layer in the electrolyte would have to be 
taken into account. 

We will use here also a description through continuum 
equations, but take care to model the actual chemistry 
and kinetics at the interface and the physical transport 
mechanisms. This results in a somewhat complicated set 
of evolution equations and we have to be satisfied with 
the more modest goal to understand whether continuum 
equations in general are able to predict the transition 
from electropolishing to pore formation. As benchmark 
for the transition we use the dispersion relation as ob- 
tained from a linear stability analysis of the flat interface 
moving at a constant velocity. If the dispersion switches 
from unstable to stable, we interpret this as the transi- 
tion from pore formation to polishing. Of course, if the 
continuum equations contain the information of the pore 
structure at all, it will not be unraveled in such a stability 
analysis. 

In Section ^ wc discuss the full nonlinear transport 
equations and their boundary conditions at the interface. 
We derive a simplified transport model which covers es- 
sential features of electrochemical etching of semiconduc- 
tors and calculate the stability of the dissolution front in 
linear order in Section III. Our results are summarized 



in Section IV. 



II. MODELING ELECTROCHEMICAL ETCHING 

A. Electrochemistry 

While anodizing silicon in hydrofluoric acid, silicon is 
dissolved in an electrochemical reaction. The detailed re- 
action mechanism is still topic of actual research. How- 
ever, during pore formation hydrogen evolution is ob- 
served whereas no hydrogen is formed during electropol- 
ishing. The valence of the chemical reaction differs in 
both cases. During pore formation the dissolution steps 
for a single silicon atom add up to Q 



Si + 6 HF + h^ 



SiF;^ 



H2+4H+ 



(1) 



For each silicon atom two elementary charges have to 
reach the interface. For electropolishing the sum reac- 
tion is HI 



Si + 6HF + 4 h^ 



SiF^ 



6H^ 



(2) 



In both cases, for pore formation and for electropol- 
ishing, the local silicon dissolution rate is proportional 
to the local electrical current density component normal 
to the interface, denoted by j^. Thus the local interface 
velocity w is 



w = 



(3) 



where F is the volume of silicon per unit area dissolved 
by a unit charge. This number is inversely proportional 
to the valence v of the chemical reaction, i.e. the num- 
ber of unit charges necessary to dissolve one silicon atom. 
The valence doubles as the current density is increased 
over the threshold value for electropolishing p^ ]. 



B. Transport equations 

The local current density is determined by the trans- 
port of reactants and reaction products in the semicon- 
ductor and the electrolyte as well as by the reaction ki- 
netics. Modeling the transport and the interface reaction 
depends strongly on the considered length scale. In n- 
doped silicon, the pore spacing is typically some microns 
whereas in p-silicon nanometer sized pores are formed. 

The shorter the considered length scale the more de- 
tailed the model has to be. The mean free path of charge 
carriers in the semiconductor is of the order of some ten 
nanometers. Transport on this length scale can be de- 
scribed by Boltzmann equations, but their nonlocality 
makes them difficult to analyze [|5| . 

On the nanometer scale, quantum effects start to play 
a role and the electrical double layer at the inter- 
face in the electrolyte, i.e. the Helmholtz layer, has to 
be taken into account |16|. At this level, details about 
the electrochemical reaction path way have to be fed into 
electronic structure calculations to determine the bound- 
ary conditions. Such detailed knowledge is not available 
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and the presence of an electrolyte makes the calculations 
even more complicated. 

To avoid these difhculties we restrict our model to 
length scales large as compared to the mean free path 
in the semiconductor, i.e. larger than 100 nm. Then, 
the transport in the semiconductor and in the electrolyte 
can be described by drift and diffusion, i.e. by Nernst- 
Planck equations. The current density of electrons, j„, 
and holes, jp, is given by 



in ^ e Dn Vn + /i„ n E, 
jp = -e ZJp Vp + p E, 



(4) 



where D^/p is the diffusion constant, ^„/p the mobility, 
e the elementary charge, and E the electric field. The 
electron and hole concentrations n and p determine the 
local charge density and thus the electric field via the 
Poisson equation 



V • E = — [p-n + N), 
esi 



(5) 



with N the density of ionized dopant atoms and egj the 
dielectric constant of silicon. For p-doped silicon N is 
negative and for n-doped silicon positive. Since the elec- 
tric field is determined by the charge carrier concentra- 
tion, the products n E and p E in Eq. represent non- 
linearities. Another nonlinearity appears due to produc- 
tion and recombination of electron hole pairs in the con- 
tinuity equation 



e pn- Paqn^ci 



-V-jp, 



(6) 



where r the life time of the charge carriers and ncq and 
Pcq is the equilibrium electron and hole concentration, 
respectively p^ . We make the quasistatic approxima- 
tion neglecting the time derivative of the concentration 
fields. This is justified if the relaxation time for any field 
involved in the dissolution process is much faster than 
the interface movement. The electrical charge is con- 
served due to V • j = V • (j„ + jp) = 0. The source term 
in Eq. (^) is derived from a law of mass action for the 
recombination reaction e^ + h+ ^ 0. 

If convection is negligible, the transport of molecules 
and ions in the electrolyte can be described analogously. 
HF, H2, H+ and SiFg~ have to be considered as well as 
fluoride F~ and 0H~ . For each component X a Nernst- 
Planck equation gives the particle current density, de- 
noted by ix (in contrast to the electric current density 



ix = -Dx VCx + qx MX Cx E, 



(7) 



with the particle's charge q^- The reactions HF ^ 
H+ + F" and H2O ^ H+ -I- 0H~ have to be taken into 
account in source terms for the continuity equations of 
the corresponding current densities 



V-i, 



OH- 



V-iH+ 



1 CoH- C'h+ — 
TH20 C^^o 
1 Cob- C'h+ - 



'''H20 



(~<2 



HF 



THF 



V • ip- = 



1 

T"HF 



CH2O 
'^-^HF 



(8) 



CH2O 

The equilibrium constants for the water dissociation and 
the HF hydration are = IQ-i^ 2^ and J^hf 
S.SlO^-^Sfi In (|) the water concentration, con- 

sidered as a natural constant of unit , has been added 
to correct the units iQ. In the electrolyte the Poisson 
equation is 

V • E = (Ch+ ~ C'p- — Cqh- — 2 C, 



SiF;^ 



(9) 



with the electrolyte's dielectric constant egio. As in the 
semiconductor, the law of mass action type source terms 
and the coupling of the ion concentrations to the electric 
field represent nonlincaritics which make the transport 
equations considerably more complicated as compared to 
equations used to describe directed solidification or vis- 
cous fingering for example |20|. Moreover, especially in 
the description of the transport in the electrolyte, there 
are poorly understood features. First of all, there are 
many other ions in the electrolyte that do not take part 
in the dissolution reaction directly but affect the electric 
field and the transport properties. Secondly, during pore 
formation, silicon enters the solution as HSiFa and re- 
acts to SiFg~ in the solution jl^. The reaction rate for 
this process is not known and it is possible that a con- 
siderable amount of HSiFa is present near the interface 
in solution. Thirdly, hydrogen bubble evolution is not 
modeled as well as convection. Nevertheless, we propose 
to proceed with the nonlinear transport equations (11) to 



C. Boundary conditions 

At the interface between the semiconductor and the 
electrolyte, the dissolution reactions Eqs. (]l|) and (|^) have 
to be taken into account. In both equations, the electric 
current density component normal to the interface (de- 
noted by the subscript _L) is given by the difference of 
the forward and backward reaction rate. In a law of 
mass action type approximation, these rates are given by 
the concentration product of the reactants and the reac- 
tion products, respectively. Thus for the pore formation 
reaction we obtain 



"Tporc (C'hF P ~ '7poro C'gjp2- C] 



H+ 



(10) 



with the reaction rate Fporc- The parameter Typoro is a 
measure for the equilibrium concentration product for 
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this reaction. The current density is taken positive for 
currents flowing from the electrolyte into the semicon- 
ductor. For reaction Eq. (||) we obtain correspondingly 

j± — — Tpolish (ChF ~ '7polisli C'siFg" ^H+) • (H) 

In both cases, the concentration of crystalline silicon is a 
constant summarized in the parameters F and rj. Surface 
tension could be accounted for by a curvature dependence 
of 77 but cannot play a role at the length scales discussed 
here. Obviously, these boundary conditions are highly 
nonlinear. The particle current density components nor- 
mal to the interface represent the stoichiometry of the 
corresponding chemical reaction. In case of Eq. @ this 
is 



l»HF± 


= 6 iH+j. , 




— Jn± , 


1 . 






= i}i2± ' 


4 . 






= *H+J_ ; 






1 




-Jp± 


= «SiF^-^ 



and for Eq. (|) 

1 , , 

-jp_L =4^sip2^^, (13) 

-Jpj_ = 4lH+_L- 

Thus the current densities of all species taking part in the 
reaction are determined by fixing the current density of 
one of them. The normal current density component of 
species which do not take part in the dissolution reaction 
is zero. 

The inner Helmholtz layer in the electrolyte is not de- 
scribed by the ion transport equations. This layer has 
a very high capacity as compared to the diffuse part of 
the electrical double layer in the electrolyte and the de- 
pletion layer in the semiconductor ||2l|. For this reason, 
the main potential drop across the interface occurs in the 
depletion layer in the semiconductor. Thus the electrical 
potential (E = — VV) is continuous at the interface 
in a first approximation. For the electric field boundary 
conditions as known from electrostatics are used, i.e. the 
tangential component of E and the normal component of 
e E are continuous across the interface. 

We do not specify the boundary conditions at the cath- 
ode or at the backside of the wafer, since no experimental 
evidence is reported for an influence of the cathode or the 
wafer backside on the dissolution process. Outside the 
depletion layer the semiconductor is electrically neutral 
and can be treated as an Ohmic conductor. The electron 
and hole concentrations have equilibrium values. In the 



electrolyte, there is more variety. Due to the consump- 
tion of reactants and the accumulation of reaction prod- 
ucts, the composition of the electrolyte changes in time. 
However, usually the electrolyte is stirred. In a large con- 
tainer, this means, that in a certain distance from the 
anode, the electrolyte is homogeneous, approximately in 
equilibrium, and has a composition which hardly changes 
in time. Thus, a reasonable boundary condition for the 
model equations is to fix the concentration of the elec- 
trolyte components at a certain distance from the inter- 
face to the equilibrium values. This distance depends 
strongly on stirring and the current density and must 
remain as a free parameter. Such kind of boundary con- 
dition is a very crude approximation, since convection 
certainly plays a role even in the diffusion layer, i.e. the 
region near the interface where the electrolyte is not ho- 
mogeneous due to the applied current. To include con- 
vection in the model is in principle possible. On the scale 
of the pores convection should not play a role because of 
the high viscosity of hydrofluoric acid. 

The model described in this section should be reason- 
ably close to the physics and chemistry of the etching 
process on length scales large compared to the mean free 
path in the semiconductor. We have to deal with nonlin- 
ear transport equations and boundary conditions. How- 
ever, the two cases, i.e. pore formation and electropol- 
ishing, are actually treated as two separate models. One 
would like to have a model that decides itself which re- 
action pathway (i.e. which law of mass action bound- 
ary condition) to take, depending on concentrations or 
current densities at the interface. Combining the two 
boundary conditions Eq. (^o|) and Eq. ([ll| ) to 

i_L = ^Tporc (C'hF P ^ ^porc C'siFg^ Ch+ ^) 

^Tpolish (C'hF 

/-'7poiishCsip2- Cg+) (14) 

gives the correct normal current density. But the stoi- 
chiometry equations (|l2|) and ( p^ ) have to be included 
too. For that one needs to know the local fraction of 
silicon atoms that are dissolved by the pore reaction (|l|) 
and the polishing reaction (^), respectively. 

Beside this constraint, analyzing the above developed 
model (analytically or numerically) would be a really de- 
manding venture. The parameter space is large and even 
though many parameters have a direct physical interpre- 
tation, experimental values are not available to fix them. 

Our goal is to study the transition from pore forma- 
tion to electropolishing, in particular to investigate the 
mechanism that can lead to such a transition. Therefore 
we study a simplified model that captures key features 
of the above description of electrochemical dissolution of 
silicon, namely the interplay of more than one field deter- 
mining the interface motion, a change of valence of the 
dissolution reaction, and law of mass action type bound- 
ary condition. 
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III. LINEAR STABILITY ANALYSIS 
A. Simplified model 

1. Model equations 

To have a tractable model and to gain some experi- 
ence we have to simplify and consider only one field, '^e, 
in the electrolyte and one field, ^'s, in the semiconduc- 
tor. These fields could be either the concentrations of 
one of the species in the electrochemical reaction, or the 
electrical potential, or a linear combination of fields as, 
e.g., the total amount of fluor per unit volume which is 
the sum Cp- + Chf- For simplicity, we will work with 
concentration fields in the following. To account for the 
interaction between the various fields we include source 
terms in the continuity equations which drive the fields 
to their equilibrium values ^''^ and ^'^^'^ in the electrolyte 
and the semiconductor, respectively. We assume only dif- 
fusive transport to keep the transport equations linear, 
which allows the stability analysis of a planar interface to 
be performed analytically. Nonlinear transport equations 
would lead to linearized equations for a flat and planar 
interface with nonconstant coefficients. The particle cur- 
rent densities in the electrolyte and the semiconductor, 
i_E and is, respectively, are assumed to be given by 



i_E = V*£; and is^-D^V^s, 



(15) 



with the diffusion constants De and Ds- The continuity 
equations are then 



V • lE/S 



-D 



E/S 



1 



E/S 



TE/S 



(*iJ/S 



^ E/S) 



(16) 



with the time constants te and Tg for the electrolyte and 
the semiconductor, respectively 

At the moving interface T = {{x^y, z)\z — h{x,y;t)} 
we assume a law of mass action type boundary condition 
for the reaction 



0, 



(17) 



similar to Eq. ( |10D or (^_l|). The setup is illustrated in 
Fig. Since we want to perform a linear stability anal- 
ysis of a flat and planar interface the restriction to sin- 
gle valued interfaces is no real limitation for our anal- 
ysis. Another choice of the interface reaction would be 
^E ^ ^'s, leading to linear law of mass type boundary 
conditions, see [Q. This type of boundary conditions 
does not lead to the desired properties. The law of mass 
action for Eq. ( [l7| ) is 



(18) 



and the stoichiometry is represented by 



At some distance dE and ds from the interface T, ei- 
ther the current density Ie/s or the field ^_e/s can be 
fixed. However, at least at one side the field has to be 
fixed to eliminate all gauge freedom. We will fix the cur- 
rent density in the semiconductor and the field in the 
electrolyte, i.e. isU=d, = I and 'i'E\z=-dE = ^d^- The 
normal velocity w of the interface is given by the normal 
current density 



-Fi 



s±\ 



(20) 



With the choice of sign in the above equation and in- 
terface reaction (|l^ ) the particle current in the semicon- 
ductor has to fiow to the interface to dissolve the silicon. 
This means that VPg has to be the hole concentration. 

A change in the valence of the dissolution reaction can 
then be modeled by a current dependent F{is±), which 
is inverse proportional to the valence cx y, i.e. the 
number of particles '^e/s needed to dissolve a certain 
amount of semiconductor material. Geometric consid- 
erations lead to the growth rate of the height function 
h{x,y;t) 









• is 



where h denotes the time derivative of h. 



(21) 



2. Linearized theory 

For a flat and planar interface hoit), the transport 
equations ( p^ become ordinary linear differential equa- 
tions with constant coefficients in the independent vari- 
able z. The solutions, which can be obtained analytically, 
depend on the boundary conditions a,t z = dE and z = ds 
and will be discussed later. They are linear combinations 
of exponentials or, in the limiting case te/s oo, affine 
functions of z. 

Now we assume a small perturbation Sh(x,y;t) of a 
certain wave length A = ^ of the interface and expand 
the fields (and correspondingly the current densities) up 
to first order in Sh, 



'fE/s = 'S%/s + S^E/s + 0{6h^), 
lE/s =-i%/s + ^^E/s + 0{Sh^). 



(22) 
(23) 



From Eq. ( pi| ) the time evolution of the perturbation can 
be derived 



5h = 



Fi4^) ( 



Sh i 



0' 



5i<- 



dlnF 



Ai 



5_L 



(*) 



(19) 



(24) 
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with the prime abbreviating the derivative with respect 
to z. The last term (*) is independent of the shape of 
the perturbation 6h. In terms of the valence ly it can be 
written as 



with uj{k) > will grow exponentially and are called un- 
stable, whereas modes with Lu{k) < are damped and 
stable. 

From the continuity equations ( |l6| ) it follows 



7" 



dhny 



di 



s± 



(25) 



D 



E/S 



S'^E/s - ^'i'E/s) = —S^E (26) 

' ^ te/s 



The growth speed of the perturbation Sh is thus pro- 
portional to Sh and Eq. (E4) can be written as Sh =^ 
Lu{k) Sh with the dispersion relation uj(k). Perturbations 



by comparing powers of Sh. The first order terms in the 
boundary conditions Eq. (^8|) and (|l^) for the perturbed 
fields S'^E/s are 



Si 



Shi%'^_ 



5h¥'E 



^7 



) i^S^s + Sh-^l' 



(27) 



(^Stsi_ + Sh 1^ = - (^StEi_ + Sh z^;) (28) 

At the lines z = (Ie and z = ds the fields satisfy Dirichlet 
boundary conditions. 



In this section we assume for simplicity that F, and thus 
the valence is independent of the current density. 



1. Infinite life time — double Laplacian growth 



B. Change of valence 

Independent of the solution of the first order equation 
(Eq), the third term (*) in the right hand side of the 



time evolution equation (24) changes sign with the cur- 
rent density |x passing through the interface if F (or 
the valence v) varies strongly enough. A change of sign 
of this term changes the sign of uj{k) for all k, making 
stable modes unstable and vice versa. A similar mech- 
anism has been proposed to explain the stability of the 
macropore front ]2^ . 

A sharp change of the valence of the electrochemical 
dissolution reaction from 2 to 4 electrons per Si atom at 
the critical current density for electropolishing has been 
found experimentally. Estimating the valence change 



■0 dlnf 



1.18, 



from flj, p. 39, Fig. 3.8] leads to cg^ 

enough to change the sign of Lu{k). 

This change of sign leads only to a stabilization of the 
interface if there have been no stable modes for lower cur- 
rent densities, since these would become unstable. More- 
over, our analysis cannot explain, why the interface re- 
mains stable for high current densities, i.e. for electropol- 
ishing, where the valence does not change any more. 



C. Transport and boundary conditions 

The linearized interface growth model as described in 
Sec. Ill A can in principle be solved analytically for all 



types of boundary conditions far away from the interface. 
Among the many possibilities we discuss some instructive 
limiting cases and infer the general behavior from them. 



The simplest case is the limit of infinite life time 
''"e/s ~^ 00. Both fields, "^^e and ^5, solve the Laplace 
equation and this simplified model is a straight forward 
extension of the well studied Laplacian growth model 
pO| . The solutions for the planar interface are then 



I De{-I + Tv) 

Ds"" V{IdE + DE-^d, 



I 
1^ 



{z + dE) + '^a 



(29) 



After solving the first order equations ( pq ) one obtains 
for the dispersion relation 



Lu{k) = 



Flk (vl,^(0)-^vl;0(0) 



(0) coth A: ds + §^ *s (0) tanh k ds 



Ds 



(30) 



For small fc, the sign of Lu(k) is basically determined by 
the ratio of the diffusivities and the current direction and 
goes quadratically to zero 



^ik)^Fl(^^^ 



1 



dEk^ + 0{k^). (31) 



In the limit of ^ cx) the dispersion relation goes lin- 
early to zero for k ^ with the same prefactor as above 
(i.e. substitute k^ by k). 

For large k, the dispersion relation saturates at 



;(fc) ~ FIT 



De^UO) 



- 1 
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Again the term ^ (o) ~ determines the sign. For 
De ^ Ds the dispersion relation is positive for all k but 
changes sign as 15b ^ Ds- In other words, the interface 
is unstable if the front propagates into the medium with 
the much lower diffusion constant. The ratio of diffusion 
constants at which the sign changes is determined by the 
other parameters, e.g., by De/ Ds = 1/50 in the example 
in Fig. 1^. The reason for the dependence of the stabil- 
ity on the diffusion constant is the following. Reactants 
from the semiconductor (e.g., holes) reach the pore tips 
first. The diffusive transport in the semiconductor thus 
destabilizes the interface. On the other hand, the reac- 
tants in the electrolyte reach the pore walls easier than 
the tips and thus stabilize the interface. If the diffusion in 
the semiconductor is much slower than in the electrolyte, 
then the growth speed is solely determined by the trans- 
port in the semiconductor and the interface is unstable. 
This is the standard DLA scenario |2^. In the opposite 
case, the transport in the electrolyte determines the local 
dissolution rate and the interface is linearly stable. This 
anti-DLA limit has been studied in ISsfl- 



The term 



_Ds*|(0)_ 
De *"(0) 



1 also has two zeros as a func- 



tion of / since the numerator is a quadratic polynomial in 
/. For De ^ Ds and a small negative / the dispersion re- 
lation is positive for all k and changes sign at the first zero 
Ii. Fig. I illustrates this property of uj{k). This means 
that the law of mass action boundary conditions provide 
a mechanism to stabilize an interface simply by changing 
the current density. The second zero I2 is smaller than 



De I'd 



where ^'^(0) changes sign. This leads to 



an unphysical pole in the dispersion relation. Since we 
interpret the fields '^e/s as concentrations, ^'^(O) has to 
be positive. 

The finite distances dE/s have basically the effect that 
they provide an infrared cutoff changing the dispersion 
relation from linearly to quadratically for /c — > (besides 
changing the numerical values of ^'^^5). Therefore we 
set (Ie 00 for the discussion of the model with finite 
life times in the electrolyte. 



2. Helmholtz equation in the electrolyte 

To study the effect of a finite life time for the dif- 
fusing species, we take te finite while keeping ts infi- 
nite. The choice is motivated by the high background 
ion concentrations in the electrolyte which provide a 
buffer/reservoir for particles. Since we set cIe — >■ 00, 
equilibrium is the only choice for the boundary condi- 
tion in the electrolyte far away from the interface, i.e. 

E ^ '^'e f*-"^ ^ ^ —00. The solution for the flat inter- 
face is then 



De 



DekeI'''^ 



De ke 



E ' 



(33) 



where ke = 1/\/De te is the reciprocal diffusion length. 
The dispersion relation for this case is 



r(4 + *^") V^iTi^-DsDE^^JV^^Tk^-nE) 

uj{k) = -FIk^ —, ^ ^ ^ -. (34) 

(^De i^s fc + r (4 + COthfc dj ^Kj, + fc2 + Ds k 



For small k, the sign of the dispersion relation is de- For ds 00, uj{k) goes linearly to zero with the same 
termined by the sign of the current density /, prefactor (i.e. substitute ds k^ by k). 

The dispersion relation has a limit for k —i- 00, which 
u}{k) ^ —F I ds k^ . (35) is a nonlinear function of the current density / 



;(fc) 



-FI 



TP + Ds ke {Ds ke 



2T^''^)I + DsTkI {Ds^t 



Dsv) 



Ds Ds KE {I + De ke ^b") 



(36) 



The sign of this limiting value changes with the current 
density in the same way as the sign of the dispersion rela- 
tion in the previous section. The difference here is, that 
the sign for small values of k is independently fixed by 



the sign of /. For / < 0, i.e. in the case of dissolution of 
the semiconductor, the dispersion relation is positive for 
small k and positive or negative for large k. The change 
in the stability of small wave length modes is illustrated 
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in Fig. ^. 

The reason why the sign of the dispersion relation for 
smaU k is given by the current density / only is, that in 
this limit the electrolyte does not influence the stability 
of the interface. The diffusion length 1/ke is a cutoff for 
the wave length up to which the electrolyte can stabi- 
lize the interface. Longer wave length perturbations are 
controlled only by the semiconductor. 



3. Finite life time in semiconductor and electrolyte 

The findings in the above section lead to the conjecture 
that the stability of long wave length perturbations is de- 
termined by the medium with the longer diffusion length. 
To verify this we assume a semi-infinite electrolyte and 
semiconductor to keep formulae simple. With this type 
of equation and boundary condition it is no longer possi- 



ble to fix the current density / but by changing the ration 
of 'I'^/^'g'^ the current through the interface can be con- 
trolled indirectly. The solutions for the flat interface are 
then 

^%{z)^j^Be^-^' + ^l\ (37) 

where B is the positive root of 

= rDsKsB^ 
+ (r (Ds Ks + De he "f"^) + Ds De US kb) B 
+ TDEnE{%^^'^^v)- (38) 

We write the dispersion relation in terms of the zeroth 
order fields and current density at the interface and we 
use the abbreviation ICe/s = \/^'e/s 



u^[k)- ^l^sA^) r{DElCE^UO) + DslCsmO))+DEDslCElCs- ^^"^^ 

For long wave length, the dispersion relation is 

^^(4'^°(0)- 4*0(0)) 
u;{k) -FTi°^{0)— 0/ ^^ k' + Oik^). (40) 



r {De he ^-^(0) + Ds US *!^(0)) + De Ds ke ks 



If the semiconductor has a much longer diffusion length 
(i.e. Ks <C ke), the leading term is stable and it is un- 
stable if the diffusion length in the electrolyte is much 
longer. The value of the ratio ke/hs at which the sign 
changes depends on the values of the other parameters. 
In the example of Fig. ^ the critical ratio is ke/^s = 1-83. 

For large k ^ oo, the dispersion relation has a limit 



u;{k) = -Fri^sAO) ( 



Ds 



De 



+ 0{-). (41) 



Like in the previous sections, the stability of small scale 
perturbations is determined by the ratio of diffusion con- 
stants. They are stable if the diffusivity in the semicon- 
ductor is much larger than in the electrolyte and unstable 
if Ds De- For the parameters in Fig. ||, the critical 
ratio of diffusion constants is Ds/De = 0.90/2.0. 

With the boundary conditions discussed in this sec- 
tion, the current through the interface is controlled by 
the ratio '^'^^/'i>°^. Like in the last section, short wave 
length perturbations can be stabilized by increasing the 
current over a certain threshold, i.e. increasing the ra- 
tio of 5'^''/*^. In the example in Fig. ^ this value is 
8.58/2.0. 

Reducing the ratio ^E'^'/^E'^ to the equilibrium value, 
i.e. for — ri, the zero order current density 1^(0) 



vanishes and the interface becomes marginally stable. 
Below that value semiconductor material is deposited and 
uj{k) changes sign for all k. In Fig. ^ the growth rate of 
the flat interface ho and the limit of the dispersion re- 
lation for large k are plotted against VE'^'^ for the same 
parameters as in Fig. showing the two sign changes of 
w(fc — > oo). The sign change at — 0.5 is accompa- 
nied by a reversion of the growth direction whereas at 
^'^'^ = 8.58 only the stability properties change. 



IV. DISCUSSION 

Motivated by the discussion of electrochemical etching 
of silicon in HF solutions we developed a simplified model 
in Sec. 



Ill 



This model has more than one field determin- 
ing the interface motion, includes a change of valence in 
the dissolution reaction, has nonlinear law of mass ac- 
tion type boundary condition, and includes a background 
reservoir for the reactants. The transport in both, the 
semiconductor and the electrolyte is diffusive. We per- 
formed a linear stability analysis of a flat interface for 
three limiting cases of a simplified model and found two 
mechanisms that can cause a stabilization of the interface 
at high current densities. 
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Firstly, a change of valence of the dissolution reaction 
with current density can stabilize the interface. The sign 
of the dispersion relation flips when the change of the 
valence with current density is large enough, namely if 
^5'_l(0) 'di^^"* \i° (0) > 1- This is true in general, in- 
dependent of the transport mechanisms in the semicon- 
ductor and electrolyte, the number of reactants and re- 
action products, and the type of boundary conditions. 
For electrochemical etching of silicon in hydrofluoric acid 
the change of valence at the transition from pore forma- 
tion at low current densities to electropolishing at high 
current densities is large enough. The stability of the 
interface at high current densities, after the valence set- 
tled at the electropolishing value, cannot be explained 
with this mechanism. However, there an oxide layer is 
formed which has to be dissolved chemically. At low cur- 
rent densities, when pores are formed, silicon is dissolved 
directly. This oxide layer reduces the diffusivity in the 
electrolyte considerably, leading to a stabilization of the 
interface. Our model shows, that an interface propagat- 
ing into a much more diffusive medium is stable (given 
that the reactants are the rate limiting species, not the 
reaction products). 

A second mechanism that at least partially stabilizes 
an interface are nonlinear law of mass action type bound- 
ary conditions. While the sign of the dispersion relation 
u}{k) at small k is determined by the ratio of diffusion 
lengths in semiconductor and electrolyte, the value at 
high fc as a function of the current density can change 
sign from positive to negative. The sign of a;(fc) for small 
k is determined by the ratio of diffusion lengths. The side 
with the larger diffusion length 1/k, i.e. the larger life 
time T, determines the stability of long wave length per- 
turbations. If the diffusion length of the semiconductor 
is much longer than the one of the electrolyte, the in- 
terface is unstable for long wavelengths. If the diffusion 
length in the electrolyte is much larger, u^k) is negative 
for small k. 

Our analysis shows, that a continuum model of sur- 
face growth, i.e. a moving boundary value problem for 
partial differential equations, can have a transition from 
linear stability to instability with increasing current den- 
sity. The two mechanisms discussed here are particularly 
interesting from a theoretical point of view in that the ef- 
fects of nonlinearities in the model equations can still be 
handled analytically through a linear stability analysis. 

It is questionable whether the transition from porous 
silicon formation to electropolishing can be described by 
one of the discussed mechanisms alone. The valence of 
the dissolution reaction of silicon does change in this 
transition. But the implementation in the model of 
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Sec. Ill A is oversimplified. A more realistic model would 
include two alternative reaction pathways for the species 
at the interface in the spirit of Eq. The analysis in 

this paper should be regarded as a stepping stone for the 
development of such models. 
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FIG. 1. Sketch of the simplified model. The semiconduc- 
tor is above the interface h{x, y) and the electrolyte below. 
Particles move to the interface from both sides and react with 
each other, dissolving the semiconductor. 



FIG. 2. Dispersion relation in case of infinite life time and 
for several values of Ds- Note the sign change at Ds = 50. 
The other parameters are set to F = 1, De — 1, ds = 10, 
ds = 20, = 1, I = -1 and = 20. 



FIG. 3. Dispersion relation in case of infinite life time and 
for various current densities. Note the sign change between 
/ = —1.8 and / — —1.95. The other parameters are F = 1, 
De = 1, Ds = 0.5, ds = 10, ds = 20, r] = 1, and = 20. 

FIG. 4. The dispersion relation in case of finite life time 
in the electrolyte and infinite life time in the semiconductor 
for various current densities. The sign of uj{k) for k oo 
changes at / = —3.38. The other parameters are F = 1, 
r = 1, De = 10, Ds = 1, ds = 10, r] = 1, K,E = 1, and 
= 1. 



FIG. 5. Dispersion relation in case of finite life time in 
semiconductor and electrolyte for different values of k_b. Note 
that the sign of u}{k) for small k changes at ke ~ 1.83 but not 
the sign for large k. The other parameters are f = 1, F = 1, 
De = 2, Ds = 1, K,s = 1, e = l, ^l/"^ = 2, and = 10. 

FIG. 6. Dispersion relation in case of finite life time in 
semiconductor and electrolyte for different values of Ds ■ Note 
that the limit of uj{k) for k —> oo changes sign at Ds = 0.90. 
The other parameters are F = 1, T — 1, De — 2, us = 1, 
HE = 10, e = 1, = 2, and = 10. 

FIG. 7. Dispersion relation in case of finite life time in 
semiconductor and electrolyte for different values of 'i'"^, i.e. 
different current densities. Note that uj{k) for k <^ 1 remains 
positive. The other parameters are F = 1, F — 1, De = 2, 
Ds = 1, Ks = 1, KE = 10, e = 1, and = 2. 



FIG. 8. The growth velocity of the ffat interface (dashed 
line) and the limiting value of Lo{k) for — + oo (dotted 
line) are plotted against 'I'^''. The growth velocity and ijj{cx)) 
change sign at = 0.5. At = 8.58 only lj(cx3) changes 
sign. The other parameters are the same as in Fig. M. 
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